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PROCEDE D'ETALONNAGE D'UN SYSTEME D'IMAGERIE, SUPPORT DE MEMOIRE ET DISPOSITIF 
ASSOCIE. 

_ Procede d'etalonnage d'un systeme d'imagerie, dans 
lequet on etalonne des parametres intrinseques relatifs aux 
differents types d'images pouvant etre prises dans des po- 
sitions angulaires predeterminees du moyen d'emission, et 
du moyen de reception, les parametres intrinseques etant 
estimes a I'aide d'un fantome bidimensionnel pour etablir un 
champ de distorsion pour une angulation particuliere, le 
champ de distorsion pour une angulation quelconque etant 
calcule par interpolation polynomiale des champs de distor- 
sion etablis chacun pour un point de controle, un point de 
controie etant defini par une angulation particuliere par axe 
de rotation du systeme d'imagerie, devant faire I'objet d'un 
etalonnage. 
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Proc6d6 d'etalonnage d f un systeme d'imagerie, support de memoire 

et dispositif associe 

La pr6sente invention concerne le domaine de l'6talonnage d'une 
chaine de traitement d'images pour corriger des distorsions, en particulier 
dans des images radiologiques bidimensionnelles, par exemple dans le 
domaine de la radiologie. 

5 De fa9on connue, les appareils de radiologie comprennent un 

moyen demission d'un faisceau de rayons X ou ensemble radiogene tel 
qu'un tube k rayons X et un moyen de reception dudit faisceau tel qu'un 
detecteur k l'etat solide ou encore un scintillateur et une camera video, par 
exemple de type CCD. 

10 Le moyen Remission et le moyen de reception de rayons X sont 

generalement supportes par un systdme mobile k un ou plusieurs axes, 
permettant des prises de vue selon diff6rentes incidences. Le moyen de 
reception du faisceau de rayons X est reli6 k des moyens de traitement 
d'images permettant la g6n6ration d'images tridimensionnelles k partir 

15 d'une serie d'images bidimensionnelles prises par le moyen de reception. 

II s'avfcre important d'obtenir une bonne correspondance entre 
une image tridimensionnelle et une image bidimensionnelle prise lors 
d'une etape au cours de laquelle le patient subit une action particuliere 
telle que la pose d'un catheter dans le domaine de l'angiographie, pour 

20 pouvoir suivre le d6placement du catheter dans l'image bidimensionnelle 
mais egalement dans l'image tridimensionnelle. 

A cet effet, un etalonnage de la machine de radiologie est 
effectu6 lors de sa mise en service puis k intervalles reguliers, pour tenir 
compte des imperfections de differents elements de la machine de 

25 radiologie, tels que la deformation des supports du moyen d'emission et du 
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moyen de reception de rayons X, les imprecisions des capteui s, les defauts 
du recepteur, etc. 

Jusqu'a pr6sent, un etelonnage pour certaines configurations de 
la machine est effectu6 a intervalles courts, par exemple une semaine, d'oii 
une immobilisation de la machine et un cout de main d'oeuvre importants. 

La presente invention propose, notamment, un procede 
d'etalonnage perfectionne permettant une precision amelioree des images 
obtenues pour toutes les configurations de distorsion et des parametres 
intrinseques de la geometrie conique de la machine. 

La presente invention propose, notamment, un proc6de 
d'etalonnage economique nSduisant la duree des operations d'etalonnage 
d'entretien effectuees regulierement au cours de la vie d une machine 
d'imagerie. 

Le procede d'etalonnage, selon un aspect de l'invention, est 
destine" a un systeme d'imagerie, comprenant un moyen d'emission d'un 
faisceau d'energie, un moyen de reception du faisceau d'energie apres 
qu'il a traverse un organe a 6tudier, un moyen de commande du moyen de 
reception, un moyen de support dudit organe et un moyen de mise en forme 
d'images. On etalonne des parametres intrinseques relatifs aux differents 
types d'images pouvant €tre prises dans une configuration predeterminee 
du systeme, les parametres intrinseques 6tant estimes a l'aide d'un 
fantome bidimensionnel pour etablir un champ de distorsion pour une 
configuration particuliere, le champ de distorsion pour une configuration 
quelconque etont calculi par combinaison d'une pluralite de champs de 
distorsion de causes supposees independantes. 

Dans un mode de realisation de l'invention, on determine un 
champ de distorsion D Q du au champ electromagnetique terrestre a partir 
d'une image dudit fantome prise dans une position centrale du systeme 
d'imagerie selon chaque axe de rotation, ledit champ de distorsion D Q 
etant independant de l'angulation. Le champ de distorsion D 0 sera 
recalcule frequemment pour tenir compte des variations du champ 
electromagnetique terrestre. Au contraire, les autres distorsions evoluent 
de fa ? on plus lente et pourront 6tre recalculees moins frequemment, par 
exemple annuellement. 

Dans un mode de realisation de l'invention, on etalonne des 
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param&tres intrinseques reiatifs aux differents types d'images pouvant 
etre prises dans des positions angulaires predetermines du moyen 
demission et du moyen de reception, les parametres intrinseques etant 
estim6s a l'aide d'un fantSme bidimensionnel pour 6tablir un champ de 
distorsion pour une angulation particuliere, le champ de distorsion pour 
une angulation quelconque 6tant calcule par interpolation polynomiale 
des champs de distorsion etablis chacun pour un point de controle, un 
point de controle etant defini par une angulation particuliere par axe de 
rotation du systfeme d'imagerie, devant faire l'objet d f un etalonnage. 

De preference, les points de controle sont d6finis de fa?on telle 
qu'ils recouvrent totalement et regulierement le domaine des orientations 
possibles de la machine consider6e. Si autour d'un axe, la rotation est 
possible sur plus de 180°, il suffira de prendre -90°, 0° et +90° comme 
angulations particulieres, 0° 6tant la position mediane par rapport aux 
extremums. 

De preference, la fonction ^interpolation est de type 
biquadratique ou bicubique. 

Avantageusement, on etablit une relation entre des parametres 
d'interpolation u, v et 1'angulation, la relation 6tant de type u(x) = ax 5 +bx 3 
+ cx, x etant une coordonnee d'un point, de fa?on que le champ de 
distorsion soit determine en fonction de 1'angulation, la fonction 
d'interpolation permettant de determiner la distorsion en fonction des 
parametres d'interpolation u, v. 

Avantageusement, la fonction u(x) est de type u(x) = oc(A+p)x 5 + 
X(A+5)x 3 + e(A+<t>)x, A etant une constante dependant du type de systeme 
d'imagerie, a, P, %> 6, e et <|) 6tant des constantes determin6es 
experimentalement 

Plus particulierement, la fonction u(x) est de type u(x) =(A- 
l,5)x 5 +(2,5-2A)x 3 + Ax. A peut etre sup6rieur ou egal k 0 et inferieur ou 
£gal h 15/8. 

De preference, on utilise un paramfetre d'interpolation pour 
chaque angle d£finissant la position du systdme d'imagerie. A sera calcule 
pour chaque angle. Les points de controle seront au nombre de (n+l) 2 , 
avec n valant au minimum 2 pour Interpolation biquadratique et 3 pour 
1'interpolation bicubique. 
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Dans un mode de realisation de l'invention, on dStermine un 
champ de distorsion Dj du a une inversion d'image par inversion du sens de 
balayage du moyen de reception du faisceau d'energie par comparaison 
d'une image dudit fantome retournee apres sa prise et d'une image dudit 
fantome prise en balayage inverse, et ce pour les changements de balayage 
utilises. 

Dans un mode de realisation de l'invention, on determine un 
champ de distorsion D R du a un changement de resolution du moyen de 
reception du faisceau d'energie par comparaison d'une image dudit 
fantdme prise avec une premiere resolution et d'une image dudit fantome 
prise avec une deuxieme resolution differente de la premiere, et ce pour 
les changements de resolution utilises. 

Dans un mode de realisation de l'invention, on determine un 
champ de distorsion D M du a un changement de matrice d'acquisition du 
moyen de reception du faisceau d'energie par comparaison d'une image 
dudit fantome prise avec une premiere matrice d'acquisition et d'une 
image dudit fantome prise avec une deuxieme matrice d'acquisition 
differente de la premiere, par exemple lors d'un changement d une matrice 
d'acquisition de 512x512 pixels a une matrice d'acquisition de 1024x1024 
pixels, et ce pour les changements de matrice d'acquisition utilises. 

L'invention propose, egalement, un support de memoire 
contenant des instructions pour mettre en ceuvre un procede d'etalonnage 
d'un systeme d'imagerie comprenant un moyen d'emission d'un faisceau 
d'energie, un moyen de reception du faisceau d'energie apres qu'il a 
traverse un organe a etudier, un moyen de commande du moyen de 
reception, un moyen de support dudit organe et un moyen de mise en forme 
d'images, procede dans lequel on etalonne des parametres intrinseques 
relatifs aux differents types d'images pouvant etre prises dans une 
configuration predetermines du systeme, les parametres intrinseques 
etant estimes a l'aide d'un fantome bidimensionnel pour etablir un champ 
de distorsion pour une configuration particuliere, le champ de distorsion 
pour une configuration quelconque etant calcule par combinaison d'une 
pluralite de champs de distorsion de causes supposees independantes. 

L'invention propose, egalement, un systeme d'imagerie 
comprenant un moyen d'emission d'un faisceau d'energie, un moyen de 
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reception du faisceau d'energie apres qu'il a traverse un organe a etudier 
un moyen de conunande du moyen de reception, et un moyen de xnise en 
forme d images. Le systeme comprend un moyen d'etalonnage des 
parametres intrinseques relatifs aux differents types d'images pouvant 
etre prises dans une configuration predetermine du systeme les 
parametres intrinseques etant estim6s a 1'aide d'un f an t6me 
b.dimensionnel pour etablir un champ de distorsion pour une 
configuration particuliere, le champ de distorsion pour une configuration 
quelconque etant calcule par combinaison d'une plura.ite de champs de 
distorsion de causes supposees independantes. 

Dans un mode de realisation de IWntion, le fantome 
mdimensionnel comprend une grille. 

Dans un mode de realisation de Invention, on effectue un 
eta,onnage lors de Installation du systeme d'imagerie. II peut S ' agir de 
1 etalonnage des parametres de la geometrie conique et de la constante A 

En outre, on peut etalonner des parametres extrinseques relatifs 
a la position dudit organe par rapport au systeme d'imagerie, du moyen de 
reception et du moyen de support, pour determiner la geometrie 
d acquisition et le champ de distorsion quelle que soit la position du 
moyen d , et du moyen de ^ ^ ^ 

etant calcules a I'aide d'un fant6me tridimensional pour des positions 
predeterminees du moyen d'emission.Le meme fantome peut servir a 
1 etalonnage des deux types de parametres. 

Avantageusement, les positions predeterminees pour le calcul 
des parametres intrinseques tiennent compte de la distance entre moyen 
d emission et moyen de reception, de la Solution, et de la macrice image 
Dans un mode de realisation de 1'invention, au moins trois 
images dun fant6me sont prises pour trois valeurs differentes de la 
distance entre moyen d'emission et moyen de reception 

Dans un mode de realisation de Invention, la position du 
systeme dimagerie est mesuree par des capteurs. 

Dans un mode de realisation de 1'invention, on effectue un 
etalonnage a intervalles de temps predetermines. 

^invention permet de corriger des distorsions dans de nombreux 
types d appareils de d'imagerie. notamment de radiologie avec des images 
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bi ou tridimensionnelles. 

Une autre application serait simplement d'obtenir ties images 
g^ometriquement correctes afin de pouvoir effectuer des mesures plus 
fiables dans les images bidimensionnelles. Ceci est egalement valable 
5 pour l'estimation des parametres intrins£ques de la g6om6trie conique qui 
permettent de deduire, entre autres, le grossissement dans les images, ce 
qui aide au recalage 2D/3D mais 6galement permet de traduire des mesures 
effectu6es dans les images (done donn£es en pixels) en des mesures 
m&riques (done en millimetres). 

10 On cherche k estimer le champ de distorsions present dans les 

images radiographiques, quelle que soit la configuration de la machine 
lors de la prise d'image. L'interet de pouvoir estimer ce champ de 
distorsion est que Ton peut gr£ce k cette connaissance corriger les images 
des distorsions. Les images radiographiques sont alors geometriquement 

15 correctes, ce qui permet d'utiliser un modele de projection conique dans 
un proc6de de recalage 2D/3D, mais encore de pouvoir faire des mesures 
plus fiables dans les images 2D. 

Un procede special est necessaire, etant donne le grand nombre 
de configurations possibles de la machine. Le principe general du proc€de 

20 est de retrouver la distorsion differentielle par rapport k un champ de 
distorsions D Q de reference. Cette distorsion differentielle est due 
principalement a quatre causes intervenant lors d'un changement de 
configuration de la machine. 

1) Changement de sens du balayage camera 

25 2) Changement de matrice d'acquisition 

3) Changement de taille de champ de vue 

4) Changement de Indentation spatiale de Tarceau portant 
T6metteur et le recepteur. 

Une configuration de la machine peut se definir comme la 
30 donn^e de 1'etat du balayage camera (par exemple : normal, inverse 
horizontalement, inverse verticalement ou inverse horizonto- 
verticalement), la taille de la matrice d'acquisition (par exemple : 256, 
512 ou 1024 pixels), la taille du champ de vue (par exemple: 11 cm, 16 cm, 
22 cm ou 30 cm) et l'orientation de Tarceau dans Tespace (donnee par deux 
35 angles, par exemple les angles anatomiques). 
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Une configuration de reference est definie en fixant : le sens du 
balayage camera (par exemple normal), la taille de la matrice 
dacquisxuon (par exemple 512 pixels), la taille du champ de vue (par 
exemple 22 cm) et 1'orientation de l arceau dans 1'espace (par exemple 
pour le plan frontal, la position anatomique correspondant a la position de 
face : angles anatomiques RAO/LAOM) et CRA/CAU=0 degres. Pour le 
Plan lateral, nous prendrons de preference la position de profil : angles 
anatomxques RAO/LAO=90 degrds et CRA/CAU=0 degres). 

^ Le champ de distorsions D Q est estim<* a partir d'une image du 
rlZZL ^ bidimenSi ° nneI aC *™ e dans configuration de 

Par conception du fantome de grille, le champ de distorsion est 
equivalent a la donnee de la distorsion en tous les noeuds, c'est-a-dire un 
vecteur bidimensionnel donnant le decalage suivant les lignes et les 
cokmnes de 1'image entre la position du noeud detecte dans 1'image de 
grille et la position theorique qu'il devrait avoir si 1'image ne contenait pas 
de distorsion. 

Un champ de distorsion differentiel est la difference noeud a 
noeud entre deux champs de distorsions et est done un champ de distorsion 
lui-rneme. 

On note »o» Implication d'un champ de distorsion a une image 
Par exemple, soit I une image et D un champ de distorsions, l image I'=D o 
I est 1 ^age I distordue par le champ de distorsion D. De meme, soit un 
autre champ D', Timage I" = D' o I' est 1'image distordue par le champ de 
distorsH>„ D Nous avons done r = D'oDoI = (D* o D) o I. Le champ de 
dorsum D • = D' o D est done le champ de distorsion qu'il faut appUquer 
k I : pour obtenir 1'image I". On dit que D' est le champ de distorsion 
differentiel entre le champ D et le champ D". 

On suppose que quatre causes de variation de distorsion 
(changement, par rapport a la configuration de reference, dans le balayage 
camera, dans la matrice d'acquisition, dans le champ de vue et/ou dans 
I orientation de 1'arceau) sont independantes. Nous pouvons alors ecrire 
tout champ de distorsion D, observe dans les images acquises pour une 
configuration particuliere de la machine comme etant compose de 
lapphcation successive du champ de reference D 0 , puis du champ 
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differentiel du au changement du balayage cam6ra D,, puis du champ 
differentiel du au changement de la matrice d'acquisition D M , puis du 
champ differentiel du au changement dans la taille du champ de vue D F 
puis du champ differentiel du au changement dans l'orientation de l'arceau 
D A , soit : 

D = D A o D F o D M o Dj o D 0 

Le procede pennet de retrouver ces champs differentiels pour 
toutes les configurations possibles de la machine. 

Concernant les trois premieres causes de variation de la 
distorsion (inversion camera, matrice d'acquisition, champ de vue), les 
valeurs que peuvent prendre ces parametres sont en tres petit nombre 
(moins de 5 sur les machines actuelles). Invoquant leur independance, on 
peut estimer le champ de distorsion differentiel du a un changement dans 
un seul d'entre eux, quelle que soit la valeur des autres parametres. En 
particulier, on peut utiliser pour les autres parametres les valeurs d6finies 
pour la configuration de reference. Par exemple, pour les changements 
dans Tinversion camera, quatre configurations d'etalonnage peuvent etre 
definies : 

. i! : tous les parametres sont les memes qu'en configuration de reference, 

sauf le balayage camera inverse horizontalement. 
. I 2 : tous les parametres sont les memes qu'en configuration de reference, 

sauf le balayage camera inverse verticalement. 
. I 3 : tous les parametres sont les memes qu'en configuration de reference, 

sauf le balayage camera inverse verticalement et horizontalement'. 
On estime les champs de distorsions D,, D 2 et D 3 en acquerant 
une image de grille dans les configurations I,, I 2 et I 3 respectivement. On 
determine ensuite les champs de distorsion differentiels D n ,D I2 etD I3 par 
comparisons de D,, D 2 et D 3 respectivement avec D Q . La comparison 
revient a etablir le vecteur de distorsion differentiel pour chaque noeud de 
la grille. 

Pour retrouver la distorsion D pour une configuration 
particuliere, il suffira de regarder l'6tat du balayage camera : 
• S'il n'y a pas d'inversion camera, la distorsion Dj sera nulle 
. Si le balayage camera se fait en inversion horizontale, D, = D n 
. Si le balayage camera se fait en inversion verticale, Dj = 
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• Si le balayage camera se fait en inversion verticale et horizontale, 
D, = D, 3 . 

On procede de maniere similaire pour etablir des champs de 
distorsion diff6rentiels correspondant aux changements de matrice 
5 d'acquisition et de taille de champ de vue. 

Le champ de distorsion, mSme lorsque l'inversion cam6ra, la 
matrice d' acquisition et le champ de vue sont gardes constants, va varier si 
Torientation de Tarceau est modifiee. Ceci est dfl au fait que la distorsion 
apparait en partie parce que les electrons qui sont cr€6s dans 
10 Tamplificateur de brillance voient leur trajectoire modifiee par Taction du 
champ electromagnetique local, cette modification etant fonction de 
Tangle entre la direction du champ electromagnetique local et la 
trajectoire des Electrons. Si Torientation de l'amplificateur de brillance 
est modifiee, Tangle respectif entre le champ electromagnetique local et 
15 la trajectoire des electrons est 6galement modifi6 et une variation de 
distorsion s'ensuit. 

A la difference des trois premiers param&tres, Torientation de 
Tarceau prend des valeurs continues (done une infinite de valeurs 
potentielles). On ne peut done pas a priori estimer un champ de distorsion 
20 diff6rentiel pour chaque orientation : on le fait pour un nombre fini de 
positions, dites de controle, et on en d6duit le champ differentiel pour 
toute autre position quand on en a besoin, et non pas pendant Tetalonnage, 
par interpolation biquadratique. 

Un mode de realisation particulier de Tinvention est illustre par 
25 les dessins annexes : 

- la figure 1 est une vue en perspective d'un appareil de 
radiologic a trois axes qui peut etre utilise pour mettre en oeuvre le 
procede ; 

- la figure 2 est une vue d'une grille d'etalonnage telle 
30 qu'observee sur un moyen de visualisation d'images associ6 au systeme de 

radiologic ; 

- la figure 3 est une courbe de la fonction u(x) pour plusieurs 
parametrages ; et 

- la figure 4 est un organigramme des 6tapes de procede. 

35 Comme on peut le voir sur la figure 1, Tappareil de radiologic 
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comprend u„ p,ed 1 en fonne de L . avec une base 2 sensiblemem 

la ba e 2. A 1 extreme opposee 5, la base 2 comprend u„ axe de rotation 
parallele au support 3 e, au.our duoue, ,e pied es. eapab.e de rounrer Un 

support 3 de facon rotative selon un axe 8. Le bras de support 6 peu, 
presenter la fortrre d'une baYonne.te. Un bras 9 circuUite en tZ. de C es 

apte a couhsser de facon rorative autour d'un axe 13 pa, rapport a 
I extremite 10 du bras de support 6. ^ 

dvw • v i 151 ^ e " ° 9 SUPPO " e tUbC 4 rayons X 11 et un detecteur 
Z 2 M T »° - '-an, ^ce. L.e 

focal du rube 1 1 au centre de la surface plane du detecteur 12 Les tro.s 

s ~ r Uon du pied du bras de supp ° n 6 - *■ - - « - 

mutuellement perpendiculaires. 

„ri,„, ' able 15 ' Pr * VUe P ° Ur recevoir ™ P*'*™. Posside une 

one„ta„o„ longitudinale alignee avec .'axe 8 en position de repos 

A trtre d'exemple, un jeu d'images d'un vaisseau sanguin e Sl 
acquts pourdifferen.es positions de ,'apparei. de radiologic. Un! image 
urd.mens.onnelle peu, etc ensuite recons,rui,e a partir du Jeu damages 
brdrmens .onnel.es. Des vues bidimensionnelles I ZZ 
trrdrmensronnelle peuven, e.re produi.es de facon in.erac.ive 

Sur la figure 2, es. momre un exemple d'une image d'une grille 
detalonnage «e.,e que qu'on peu, .'observer sur un ecran de visuahsSln 
d un apparel, de radiologic. La grille es. disposee a la place d'un obt, a 

uZT'ir 7 le r jet du faisceau * ™>™ * — >< — » x 

11 et le detecteur de rayons X 19 t * * 

Tt , Q ,. . yons A 12 - La grille est a mai les carrees en 

P ~ r::;:. relaMvem * m opa - - — * - — - 

On consta.e que I'image de la grille n'en constitue pas la 

ZZ ^ U d Ct,0nf,d41eSi,eC — el ^^-conve„ab,e m e„, pluU 
les bords presen.en. des dis.orsions que ,'on cherche a calcu.er En un 
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point donne de la grille, la distorsion peut etre caract6ris6e par un vecteur 
ayant une direction et un module. 

Parmi les causes de distorsion, on constate que certaines, comme 
les d6fauts du detecteur, sont indSpendantes des angulations autour des 

5 axes de Tappareil de radiologic tandis que d'autres dependent desdites 
angulations. Un etalonnage id£al n6cessiterait de calculer la distorsion en 
tout point du champ image, pour toute angulation selon chacun des axes, et 
plus generalement pour toute configuration de la machine ce qui 
n^cessiterait d'6normes ressources de calcul. 

10 On cherche & obtenir un etalonnage satisfaisant a cout 

raisonnable et relativement rapidement. Dans ce but, on effectue un calcul 
du champ de distorsion pour quelques angulations dites de controle et on 
en d6duit un champ de distorsion pour les autres angulations. Pour 
simplifier les calculs, on considSrera la position spatiale du tube 11 et du 

15 detecteur 12, non par rapport aux axes de rotation de 1'appareil de 
radiologic, mais par rapport & deux axes, de preference orthogonaux entre 
eux et passant par l'isocentre, perpendiculaires & la direction du faisceau 
de rayons X, en d'autres termes definissant un plan paraltele a l'axe 13. Le 
troisieme axe d'un repere forme a partir desdits deux axes est la direction 

20 du faisceau de rayons X. Une rotation du tube 1 1 et du detecteur 1 2 autour 
de ce troisieme axe n'a pas d'effet sur la distorsion. Les angulations par 
rapport auxdits deux axes sont notees 0 et cj). 

La prise d'images peut se modeliser geometriquement comme 
une projection conique perturb6e par un champ de distorsion 

25 bidimensionnel dans l'image. 

Le detecteur d'image 12 est sensible au champ 
electromagnetique environnant, ce qui introduit des distorsions de 
l'image par rapport au modele de projection conique. Une fonction de 
projection est done vue comme la combinaison d'une projection conique et 

30 d'une fonction de distorsion qui vont dependre toutes deux de Tetat 
courant de 1'appareil. La distorsion est une fonction du plan image dans 
lui-m€me. On suppose qu'elle est suffisamment reguliere. On peut alors 
l'estirner de fa9on precise en la calculant sur un nombre de points 6pars 
dans l'image, sa valeur sur les autres points en etant deduite par 

35 interpolation. 
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Le fantdme utilise pour l'etalonnage des distorsions consiste en 
une grille. Chaque intersection de la grille est appelee noeud. La position 
de chaque noeud est detectee dans 1'image. La distorsion a chaque noeud 
est definie comme l'ecart entre sa position dans l'image et celle qu'il doit 
avo ir dans le modele ideal de la grille. Quand la distorsion est corrigee 
1'hypothese de projection conique est applicable. Cette derniere est une 
fonction homographique de l'espace tridimensional dans lequel est 
expnmS l'objet vers le plan image. Elle se traduit par 1'estimation d'une 
matnce 3 x 4 a une constante multiplicative pres. Nous avons ainsi onze 
parametres a estimer, cinq parametres intrinseques modelisam la chaine 
image et six parametres extrinseques donnant la position respective de 
l'ensemble emetteur-recepteur de rayons X et de l'objet image dans 
1 espace. Cette estimation de la projection conique est effectu6e grace a un 
autre fantome consistant en un cylindre en plastique transparent aux 
rayons X sur lequel sont portees des billes en plomb disposees en helice 
Cette helice permet d'avoir une trentaine de points visibles dans 1'image 
La detect,on de ces points et leur appariement avec ceux du modele parfait 
de 1'hehce permet de determiner les parametres de la projection conique. 

En connaissant a priori la configuration de la machine dans 
Iaquelle sera prise une image, on peut alors determiner avec une grande 
pr<5c,sion la fonction de distorsion grace a la grille et la projection conique 
grace a l'h61ice. Par approximation, on considere que le champ 
magnetique terrestre est localement uniforme et varie peu dans un 
intervalle de temps de l'ordre de quelques jours. On considere egalement 
que le champ magn6tique terrestre est regulier, ce qui permet de 
considerer un faible nombre d'echantillons ou positions de controle, done 
un faible nombre d'images de la grille. 

Pour chaque axe definissant Torientation, trois angulations de 
controle, voire quatre, s'averent utiles. II en resulte done neuf, voire seize 
positions de controle definies par une angulation et une interpolation de 
type biquadratique, voire bicubique. L'angulation consideree peut etre 
une angulation "machine" par rapport a chacun des deux axes ou 
preferablement, une angulation anatomique avec deux angles (9 <|>) 
differents des angles "machine" mais pouvant s'exprimer en fonction 
d'eux. 
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Les positions de controle seront choisies regulierement 
espacees et sym6triques par rapport a l'origine, afin d'optimiser la densite 
des positions dans Tespace parametre. Dans le cas de Interpolation 
biquadratique, un ensemble de positions sera done donn£ par des valeurs 
positives 0 O et <{> 0 , qui definissent des points de controle dans l'espace (6, 

40 : (0,0), (-e 0 , o), (6 0 , o), (o, -<» 0 ), (o, c|> 0 ), (-e 0 , ^ 0 ), (-e 0 , (e 0 , -(t> 0 ), et 

(Oc ^o)- Un nombre plus eleve de points de contrSle peut etre prevu. 

Pour obtenir une precision satisfaisante, il convient d'eviter 
Textrapolation, e'est-a-dire la generation de valeur de distorsion pour des 
angulations situees en-dehors de l'intervalle compris entre - 8 0 et 9 0 , 
- <t>o et <t>o- ° n choisira done les valeurs de 0 O et de $ Q correspondant k la 
totalite ou h une large plage de la gamme d'angulations possible de fagon 
que l'ensemble des estimations soit effectue par interpolation, 9 0 et (J> 0 
restant infdrieurs ou egaux a 90°. 

Le principe de Tinterpolation quadratique est le suivant. Trois 
points dits de controle P { = y.) dans le plan peuvent etre relies par une 
courbe quadratique : (P Q ) 

P(u) = (u 2 , u, UMCPj) 

(P 2 ) 

avec M une matrice caract&istique de Tinterpolation quadratique : 
2-4 2 
M = -3 4-1 
1 0 1 

et u le parametre de Tinterpolation variant entre 0 et 1. 

Pour Tinterpolation biquadratique, une trajectoire d'un point 
dans le plan est param6tr6e par deux parametres ^interpolation u et v, et 
s'exprime : 

( p O.O P 0,lP 0 ,2) (u 2 ) 
P(u, v) = (v 2 , v, l)M(P 10 Pj A P l 2 )M'(u) 

La generalisation a un nombre de points plus 6lev6 est possible. 
La fonction d'interpolation relie les points de controle k tout point de la 
trajectoire par les parametres d'interpolation u et v. Plus precisement nous 
relions ces points a des angles anatomiques definis en fonction d'angles de 
la machine. On peut prendre u et v comme coordonnees reduites des angles 
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anatomiques, du genre u = (© _ 9 m in)^max~®min)- Toutefois, on pr6f£rera 
prendre u = f((6^ m J/(6 max -e min )) et v = g(«H^< WW) 

f et g sont des fonctions impaires, croissantes et telle que f(l) = 1 
et f (1) ss 0, ce qui aboutit & une fonction polynomiale de type ax 5 + bx 3 + 
5 ex. 

Plus pr£cis6ment, on parvient h utiliser une fonction incluant 
une constante A qui depend du type d'appareil de radiologic et qui sera 
conservee pour une meme serie d'appareils. La fonction de distorsion 
s'£crit alors : type f(x) = <x(A+P)x 5 + %(A+8)x 3 + e(A+({>)x. 

10 Une bonne approximation est fournie par la fonction suivante : 

f(x) = (A - l,5)x 5 + (2,5 - 2A) x 3 + Ax, avec A compris entre 0 et 15/8. La 
valeur de A = 1,25 est particuli&rement int6ressante dans la mesure ou elle 
permet d'obtenir une estimation convenable de la distorsion. A titre 
d'exemple, pour un champ de 16, l'erreur moyenne est d'environ 0,8 pixel 

15 avec 1,95 pixels d' erreur maximale. Pour un champ de 22, l'erreur 

moyenne est d'environ 0,75 pixel avec une erreur maximale de 1 ,65 pixels. 
Pour un champ de 30, Terreur moyenne est de l'ordre de 0,6 pixel avec une 
erreur maximale de 1,42 pixels. 

Plus generalement, on definit N points de controle, d'ou N 

20 orientations (A^^n, chaque orientation A n etant definie par les angles 

(9 n ,<}) n ) coordonnees du point de controle numero n. N configurations de la 
machine sont d6finies : 

- le balayage camera est le m8me que pour la configuration de reference 
(par exemple balayage normal) 

25 - la matrice d'acquisition est la mSme que pour la configuration de 
reference (par exemple matrice 512) 

- la taille du champ de vue est la meme que pour la configuration de 
reference (par exemple champ de 22 cm) 

- Forientation est celle correspondant au point de controle : A n = (©n*^ 
30 On acquiert une image de grille pour ces N configurations et on 

en d6duit N champs de distorsions differentiels (D An ) ] < n < N par 
comparaisons des champs de distorsion observes pour les N 
configurations et le champ de distorsion pour la configuration de 
r6f6rence. 

35 Par interpolation, on retrouve la distorsion diff6rentielle pour 
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une configuration identique & la configuration de reference, si ce n'est que 
Indentation est differente. Soit D A la distorsion differentielle pour cette 
configuration. On determine le rectangle forme par les 9 points de 
controle, pour Interpolation biquadratique, 16 dans le cas de 
Interpolation bicubique qui entourent le point de coordonn6es (6,40 dans 
Tespace des orientations. Ces 9 (ou 16) points ne sont pas uniques en 
g6n6ral, mais les propri6t6s de continuite des interpolations choisies 
garantissent que le resultat ne d6pendra pas de l'ensemble de points 
retenus, d£s lors qu'ils encadrent bien le point (9,<I)). Nous retenons done 9 
(ou 16) champs de distorsion differentiels. 

Un champ de distorsion difterentiel est un ensemble de vecteurs 
bidimensionnels donnant la distorsion pour tous les noeuds de la grille. 

Par interpolation, connaissant rorientation courante (0,<j)), on 
retrouve pour chaque noeud de la grille la valeur de la distorsion 
differentielle d. Ainsi, on determine le champ de distorsion differentiel 
D A pour Torientation courante (6,<|>). 

Sur la figure 3, on voit les differentes allures que prend la 
fonction f(x, A) en faisant varier A entre 0 et 1,875. 

Pour parfaire l'estimation de la distorsion, on peut considerer 
une fonction d'interpolation differente selon chaque axe, par exemple en 
adaptant la constante A. La constante A peut faire l'objet d'une estimation 
automatique en fonction des donnees d'un etalonnage composees de neuf 
images. 

Dans un mode de realisation particulier de l'invention, on 
cherche a estimer deux champs de distorsion, le champ D 0 lie au champ 
61ectromagn6tique terrestre a 1'endroit ou se trouve l'appareil de 
radiologic et le champ D v H6 aux differentes angulations de l'appareil de 
radiologic On pourra ainsi effectuer un etalonnage hebdomadaire ou 
mensuel du champ de distorsion D Q en raison de Involution relativement 
rapide du champ 61ectromagn6tique terrestre et n'effectuer un etalonnage 
du champ de distorsion D v qu'h intervalles plus larges, par exemple 
semestriel ou annuel en raison de la derive plus lente de cette distorsion. 

On peut 6galement calculer un champ de distorsion d'inversion 
Dj qui est lie a une inversion de Timage par inversion du balay age video du 
detecteur. En effet, le detecteur n'etant pas parfait, Tin version du 
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balayage cree aussi une distorsion qu'il est utile d'apprehender 
notamment dans des applications de neuroradiologie, en particulier avec 
des machines de radiologic biplan equipees de deux tubes ct de deux 
detecteurs pour lesquels une inversion mecanique de la prise de vue en 
passant, par exemple, le tube de la droite a la gauche et le detecteur de la 
gauche a la droite, n'est pas possible. L'estimation du champ de distorsion 
Dj est effectuee en comparant une image obtenue en balayage normal puis 
retournee et une image obtenue en balayage inverse de la meme grille. 

En outre, on peut ameliorer l'estimation de la distorsion en 
tenant compte du changement de champ de 1'image. On estime alors un 
champ de distorsion D R lie au reglage de resolution du detecteur. 

Ainsi, un etalonnage complet est mis en ceuvre lors de la 
r6ception et de la mise en service de l'appareil. Puis un etalonnage partiel 
est mis en ceuvre a intervalles temporels courts et un etalonnage complet 
est mis en ceuvre a intervalles temporels longs. L'etalonnage partiel 
pourra consister a seulement recalculer le champ de distorsion D 0 lie au 
champ electromagnetique terrestre. 

L'etalonnage complet pourra comprendre les etapes suivantes 
voir figure 4 : 

etape 16 : installation de la grille sur le trajet du faisceau 
denergie par un responsable de maintenance, 

etape 17 : prise de neuf images de la grille pour les points de 

rr e ?ra (0 ;?' ( ~ G °' 0) ' (9 °' 0) ' ( °' ^ (0 ' *° X <^ ^ *6>. Ce* - 

%)> et (9 0 , (}) 0 ). 

etape 18 : calcul des champs de distorsion correspondants aux 
points de controle, et 

etape 19 : calcul des champs de distorsion D 0 et D v . 

II convient de noter que les etapes 17 et 18 peuvent etre 
effectu6es partiellement en meme temps pour reduire la duree de 
1 etalonnage complet. Les autres champs de distorsion D, et D R feront le 
cas echeant, 1'objet d'6tapes specifiques, telles que pour le champ D, prise 
d une image de la grille pour le point de controle (0,0) en balayage inverse 
retournement de 1'image de la grille pour le point de controle (0,0) en 
balayage normal et comparison, et pour le champ D R prise d images de la 
grille pour le point de controle (0,0) sous differentes resolutions et 
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comparaison. 

L'etalonnage partiel pourra comprendre les etapes saiv antes : 

a) installation de la grille sur le trajet du faisceau d'energie par 
5 un responsable de maintenance, 

b) prise d'une image de la grille pour le point de controle (0,0). 

c) calcul du champ de distorsion D 0 correspondant. 

On peut ensuite effectuer une combinaison du champ de 
distorsion D Q etdes champs de distorsion D v , D £> D R , et D F non modifies 
10 lors de T6talonnage partiel pour obtenir les champs de distorsion 
interpoles pour toute angulation (0, <()) lorsque Ton desire corriger une 
image de la distorsion qu'elle contient. 

Reconsiderons & present le problfeme initial qui consiste h 
estimer le champ de distorsion D pour une configuration de la machine 
15 donnee. 

On connait le champ de distorsion differentiel Dj correspondant au 
changement de balayage camera entre la configuration de reference et la 
configuration 6tudiee ("configuration courante") 

On connait le champ de distorsion differentiel D M correspondant au 
20 changement de matrice d'acquisition entre la configuration de reference et 

la configuration courante 

On connait le champ de distorsion differentiel Dp correspondant au 
changement de taille du champ de vue entre la configuration de reference 
et la configuration courante 
25 . On connait le champ de distorsion diff6rentiel D A correspondant au 

changement d'orientation de l'arceau entre la configuration de reference et 
la configuration courante (par interpolation). 

Ainsi, on peut en deduire le champ de distorsion pour la 
configuration courante de la machine : 
30 D = D A o D F o D M o Dj o D 0 

Grace a Tinvention, on dispose d'un 6talonnage economique, 
avec des durees d'immobilisation de 1'appareiJ d'imagerie en vue de la 
maintenance r^duite, et de haute precision dont une partie est realisee lors 
de Installation de l'appareil d'imagerie et n'a plus h Stre mise en oeuvre 
35 par la suite et qui est particulierement adaptee pour la mise en 
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correspondance d'images bidimensionnelles et d'images 
tridimensionnelles, qui permet les reconstructions tridimensionnelles a 
partir d'un nombre limits de vues. On peut supprimer Tutilisation de 
marqueurs externes, en radiologie, au profit d'un etalonnage effectue 
avant l'arrivee du patient, ce qui reduit la dur6e de l'examen radiologique 
effectu6 sur le patient. On peut effectuer des mesures dans les images 
bidimensionnelles avec une precision accrue. 
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REVENDICATIONS 

1. Procede d'etalonnage d'un systeme de d'imagerie, comprenant 
un moyen demission d'un faisceau d*6nergie, un moyen de reception du 
faisceau d'6nergie apres qu'il a traverse un organe k etudier, un moyen de 
commande du moyen de reception, un moyen de support dudit organe et un 
moyen de mise en forme d'images, dans lequel on etalonne des parametres 
intrinseques relatifs aux differents types d'images pouvant etre prises 
dans une configuration predetermine du systeme, les parametres 
intrinsdques etant estimes a l'aide d'un fantome bidimensionnel pour 
etablir un champ de distorsion pour une configuration particuliere, le 
champ de distorsion pour une configuration quelconque etant calcule par 
combinaison d'une plurality de champs de distorsion de causes supposes 
independantes. 

2. Procede selon la revendication 1, dans lequel dans lequel on 
determine un champ de distorsion D Q du au champ electromagnetique 
terrestre a partir d'une image dudit fantome prise dans une position 
centraJe du systeme d'imagerie selon chaque axe de rotation, ledit champ 
de distorsion D Q <Stant independant de l'angulation. 

3. Procede selon la revendication 1 ou 2, dans lequel on etalonne 
des parametres intrinseques relatifs aux differents types d'images pouvant 
etre prises dans des positions angulaires predetermines du moyen 
d'emission et du moyen de reception, les parametres intrinseques 6tant 
estimes a l'aide d'un fantdme bidimensionnel pour etablir un champ de 
distorsion pour une angulation particuliere, le champ de distorsion pour 
une angulation quelconque etant calcule par interpolation polynomial 
des champs de distorsion etablis chacun pour un point de controle, un 
point de contrSle 6tant defini par une angulation particuliere par axe de 
rotation du systeme d'imagerie, devant faire 1'objet d'un etalonnage. 

4. Procede selon la revendication 3, dans lequel un point de 
controle est proche de chaque extremum de position angulaire. 

5. Procede selon la revendication 3 ou 4, dans lequel la fonction 
d'interpolation est de type biquadratique ou bicubique. 

6. Proced6 selon l'une quelconque des revendications 
prec6dentes, dans lequel on 6tablit une relation entre des parametres 
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d'interpolation u, v et l'angulation, la relation 6tant de type u(x) = ax 5 +bx 3 
+ cx, x dtant une coordonnee d'un point, de fafon que le champ de 
distorsion soit determine en fonction de l'angulation, la fonction 
d'interpolation permettant de determiner la distorsion en fonction des 
5 parametres d'interpolation u, v. 

7. Proced6 selon la revendication 6, dans lequel la fonction u(x) 
est de type u(x) = a(A+p)x 5 -f x(A+S)x 3 + e(A+<|>)x, A 6tant une constante 
dependant du type de syst&me d'imagerie, a, (3, 5, e et (J> etant des 
constantes determin6es exp£rimentalement. 

10 8. Proc6d6 selon Tune quelconque des revendications 

pr^cedentes, dans lequel on determine un champ de distorsion Dj du a une 
inversion d'image par inversion du sens de balayage du moyen de 
reception du faisceau d'6nergie par comparaison d'une image dudit 
. fant6me retournee apr&s sa prise et d'une image dudit fantome prise en 

15 balayage inverse, et ce pour les changements de balayage utilis6s. 

9. Proc6de selon Tune quelconque des revendications 
pr6cedentes, dans lequel on determine un champ de distorsion D R du a un 
changement de resolution du moyen de reception du faisceau d'energie par 
comparaison d'une image dudit fantome prise avec une premiere 

20 resolution et d'une image dudit fantome prise avec une deuxieme 
resolution diff£rente de la premiere, et ce pour les changements de 
resolution utilises. 

10. Proced6 selon Tune quelconque des revendications 
precedentes, dans lequel on determine un champ de distorsion D M dti k un 

25 changement de matrice d'acquisition du moyen de reception du faisceau 
d'energie par comparaison d'une image dudit fantome prise avec une 
premiere matrice d'acquisition et d'une image dudit fantome prise avec 
une deuxieme matrice d'acquisition diffSrente de la premiere, et ce pour 
les changements de matrice d'acquisition utilises. 

30 11. Support de memoire contenant des instructions pour mettre 

en ceu vre un proc£d£ d'etalonnage d'un systeme d'imagerie comprenant un 
moyen demission d'un faisceau d'energie, un moyen de reception du 
faisceau d'energie apres qu'il a traverse un organe k etudier, un moyen de 
commande du moyen de reception, un moyen de support dudit org£ne et un 

35 moyen de mise en forme d'images, proc6d£ dans lequel on etalonne des 
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param&tres intrins&ques relatifs aux differents types damages pouvant 
etre prises dans une configuration predeterminee du syst£me, les 
parametres intrinseques 6tant estim6s k l'aide d'un fantSme 
bidimensionnel pour 6tablir un champ de distorsion pour une 

5 configuration particuli&re, le champ de distorsion pour une configuration 
quelconque etant calculi par combinaison d'une plurality de champs de 
distorsion de causes supposees ind6pendantes. 

12. Sy steme d'imagerie, comprenant un moyen d'6mission (11) 
d'un faisceau d'energie, un moyen de r6ception (12) du faisceau d'energie 

10 apres qu'il a traverse un organe a etudier, un moyen de commande du 
moyen de reception, et un moyen de mise en forme d'images, caracterise 
par le fait qu'il comprend un moyen d'etalonnage des paramfetres 
intrinseques relatifs aux differents types d'images pouvant etre prises 
dans une configuration pr£d6terminee du systeme, les parametres 

15 intrinseques 6tant estimes k l'aide d'un fantome bidimensionnel pour 
6tablir un champ de distorsion pour une configuration particuliere, le 
champ de distorsion pour une configuration quelconque etant calcule par 
combinaison d'une pluralite de champs de distorsion de causes suppos6es 
independantes. 
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